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Abstract 

We have examined trapping of two-armed (m = 2) nearly vertical oscillations 
(vertical p-mode) in vertically isothermal (c s = const.) relativistic disks with toroidal 
magnetic fields. The magnetic fields are stratified so that the Alfven speed, ca, is 
constant in the vertical direction. The ratio of c\/c 2 in the vertical direction is taken 
as a parameter examining the effects of magnetic fields on wave trapping. We find 
that the two-armed nearly vertical oscillations are trapped in the inner region of 
disks and their frequencies decrease with increase of c\/c 2 . The trapped regions of 
the fundamental (n = 1) and the first-overtone (n = 2) are narrow (less than the length 
of the Schwarzschild radius, r g ) and their frequencies are relatively high (on the order 
of the angular frequency of disk rotation in the inner region). On contrast to this, 
the second-overtone (n = 3) are trapped in a wide region (a few times r g ), and their 
frequencies are low and tend to zero in the limit of c 2 A /c 2 = 2.0. 

Key words: accretion, accrection disks — quasi-periodic oscillations — neutron 
stars — two-armed disk oscillations — X-rays; stars 



1. Introduction 

Discoseismology is one of important fields in studying astrophysical disks, since in some 
of them, quasi-periodic oscillations (QPOs) have been observed and most of them seem to 
be attributed to disk oscillations. In low mass X-ray binaries where the central sources are 
neutron stars or black holes, for example, QPOs are often observed. They are classified into 
high-frequency QPOs (> 100 Hz), low-frequency QPOs and low-frequency complex (0.01 ~ 
100 Hz), power-law components, and others (van der Klis 2004). The high-frequency and low- 
frequency QPOs will be, more or less, related to oscillation phenomena in the inner region of 
relativistic disks. 

A first step to examine whether a disk oscillation mode can describe the observed QPOs is 
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to compare the frequencies resulting from the mode with those of observed QPOs. In disks which 
extend far outside, discrete frequencies of oscillations will be expected only when oscillatory 
pertubations are trapped in a particular finite region of disks. In this sense, examination of 
trapping of disk oscillations and of their frequencies is an important subject in discoseismology. 

In geometrically thin relativistic disks, there are many kinds of disk oscillation modes. 
Here, we classify them into four classes in terms of node numbers in the vertical direction and 
the frequencies in the corotating frame, i.e., p-, g-, c-, and vertical p-mode oscillations (for 
details, see Kato et al. 2008; Kato 2001). (i) The p-mode is an inertial-acoustic mode with no 
node in the vertical direction 1 (nearly horizotal oscillations). The square of their frequencies in 
corotationg frame, (u — mQ) 2 , is larger than the square of the horizontal epicyclic frequency, 
n 2 , i.e., (u — mQ) 2 > n 2 , where u is the frequency of oscillation in the inertial frame, m is 
the azimuthal wavenumber and Q is the angular velocity of disk rotation, (ii) The g-mode 
oscillation (inertial mode or r-mode) has at least one node in the vetical direction (n > 1) and 
(u — mQ) 2 is smaller than k 2 , i.e., (u — mQ) 2 < hi 2 , (iii) The vertical p-mode oscillations are 
those that have at least one node in the vertical direcin (n > 1) and have high frequencies in the 
corotating frame in the sense that (u — mQ) 2 > Q 2 ± , where Q± is the vertical epicyclic frequency 
and always larger than k. This mode corresponds to "breathing mode" of Blaes et al. (2006) 
in oscillations of tori, (iv) Among oscillations formally belonging to (u — mQ) 2 > Q\, the one- 
armed (m = 1) one with one node (n = 1) in the vertical direction has a particular position. The 
oscillation is called c-mode (corrugation mode). It is nearly incompressible motions changing 
the disk plane up and down with a corrugation pattern, corresponding to warp (or tilt). 

Trapping of p-, g-, and c-mode oscillations have been extensively studied, e.g., by Kato 
and Fukue (1980) and Ortega-Rodriguez et al. (2002) for p-mode oscillations; by Okazaki et al. 
(1987), Nowak and Wagoner (1992) and Perez et al. (1997) for g-mode oscillations; and by Kato 
(1990) and Silbergleit et al. (2001) for c-mode oscillations. For reviews, see Wagoner (1999) and 
Kato (2001). Recently, much development has been made in examination of effects of corotation 
resonance and magnetic fields on trapping and excitation (or damping) of oscillations. That is, 
amplification of non-axisymmetric p-mode oscillations by corotation resonance was found by 
Lai and Tsang (2009) and Tsang and Lai (2009b), which will be a refinding of the Papaloizou- 
Pringle instability in different situations. Different from the p-mode, non-axisymmetric g-mode 
and c-mode are havily damped by corotation resonance (Kato 2003, Li et al. 2003, Latter 
and Balbus 2009 for g-mode, and Tsang and Lai 2009b for c-mode). Magnetic fields also have 
non-negligible effects on trapping. Fu and Lai (2009) showed that magnetic fields act so as to 
destroy the general relativistic self-trapping of axisymmetric g-mode oscillations. Considering 
these recent developments, Lai and Tsang (2009) suggested that among discoseismic modes a 

1 When we mention the node number in the vertical direction, it is the number of node of density perturbation 
in the vertical direction. The node number associated with vertical component of velocity is smaller than 
that of density perturbation by one. 
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possible candidate of high-frequency QPOs will be axisymmetric p-mode oscillations. 

Compared with the above many studies on p-, g- and c-mode oscillations, trapping of 
vertical p-mode oscillations are little examined. Examination of trapped vertical p-mode oscil- 
lations is, however, interesting in relation to QPOs, since their frequencies are in a wide range 
by differences of i) node numbers in the vertical and radial directions, and ii) disk parameters. 
Furthermore, the general relativity is not essential in the trapping of vertical p-mode oscilla- 
tions. Hence, the nearly vertical oscillations will be one of good candidates to describe the 
QPOs in various disks (disks of low-mass X-ray binaries to those of dwarf-navae) by a unified 
model. 

It should be noted that the interesting oscillations of vertical p-mode will be only those 
of two-armed ones, i.e., m — 2. This is because there will be no trapping in oscillations of 
m = l, and because for oscillations of m > 3 the frequencies of trapped oscillations are too high 
except for cases of large n. The oscillations with large n, however, are not interesting from 
the observational points of view. Considering these situations we restrict our attention only to 
two-armed (m = 2) oscillations with n — 1 to n — 3. 

In previous paper (Kato 2010, paper I), we have examined trapping of two-armed ver- 
tical p-mode oscillations in disks with polytopic gas, and showed how frequencies of trapped 
oscillations depend on a change of polytropic index. In this paper we restrict our attention 
to vertically isothermal disks for simplicity, but assume that the disks are subject to toroidal 
magnetic fields and examine how the frequency of trapped oscillations and the trapped region 
depend on the strength of magnetic fields. 

2. Unperturbed Disks and Equations Describing Disk Oscillations 

We consider geometrically thin, relativistic disks. For mathematical simplicity, however, 
the effects of general relativity are taken into account only when we cosider radial distributions 
of f2(r), re(r), and Q±(r), which are, in turn, the angular velocity of disk rotation, the epicyclic 
frequencies in the radial and vertical directions. Except for them, the Newtonial formulations 
are adopted. Since geometrically thin disks are considered, Q, is approximated to be the 
relativistic Keplerian angular velocity, Q^(r), when its numerical values are necessary. Here, r 
is the radial coordinate of cylindrical ones (r,(p,z), where the z-axis is perpendicular to the disk 
plane and its origin is the disk center. Functional forms of f^O"), K { r ), and are given in 

many literatures (e.g., Kato et al. 2008). 

2.1. Unperturbed Disks with Toroidal Magnetic Fields 

The equilibrium disks are axisymmetric with toroidal magnetic fields. The fields are 
assumed to be purely toroidal with no poloidal component: 




(1) 
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We further assume that the gas is isothermal in the vertical direction and the magnetic fields -Bo 
are distributed in such a way that the Alfven speed ca is constant in the vertical direction, i.e., 



(B6747rpc 
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const, in the vertical direction, where po(r,z) is the density in the unperturbed 



disks. Furthermore, the rotation is assumed to be cylindrical; e.g., the angular velocity of disk 
rotation, Q, is only a function of r. 

We now consider the vertical structure of the disks. The hydrostatic balance in the 
vertical direction is given by 
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Poc 2 , Cg(r) being the isothermal 



where p (r,z) is the pressure and related to p Q (r,z) by p 
acoustic speed. 

Since both c s and ca are constant in the vertical direction, equation (2) can be integrated 

to give 



2H 2 



po{r,z) = p 00 {r)exp 

where the scale height H is related to c s , ca, and Q± by 
c 2 + 4/2 



and B (r,z) = B 00 (r)exp 
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2.2. Equations Describing Disk Oscillations 

Now, small-amplitude perturbations are superposed on the equilibrium disk described 
above. The velocity perturbation over rotation is denoted by (u r , u^, u z ), and the perturbed 
part of magnetic field over the unperturbed one by (b r , b^, b z ). Then, the r-, tp-, and z- 
components of the linearized equation of motions are written, respectively, as 
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where p\ and p\ denote the perturbed parts of the pressure and density, respectively. Similarly, 
the induction equation gives 
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The equation of continuity is 

Another relation that we need here is a relation between p 1 and p\. Considering isothermal 
perturbations, we adopt 

Pi = Picl (12) 

Here, the azimuthal and time dependences of the perturbed quantities are taken to be 
proportional to exp[i(ojt — rmp)], where u and m are the frequency and the azimuthal wavenum- 
ber of the perturbations, respectively. The perturbations are assumed to be local in the sense 
that their characteristic radial wavelength, A, is shorter than the characteristic radial scale of 
disks, Ad, i.e., A < A D . By using this approximation, we neglect such quantities as dlnp 00 /d\nr, 
dlnBQ j dlnr , dlnH/dlnr, and <ilnf2/<ilnr, compared with terms of the order of r/A. Then, the 
r-, </?-, and ^-components of equation of motion, equations (5) - (7), are reduced to 

i (u-mn)u r -2Qu^ = -^-cl—^J, (13) 

K 2 

i(uj — mQ)u v + — u r = 0, (14) 

-^-(l + 4>-*(s-S')(^- 1 7*(5)- (15) 

In the above equations, hi defined by hi =pi/po = c 2 pi/po has been introduced. Similarly, the 
r-, <f-, and z-components of induction equation, equations (8) - (10), are reduced to 

i{u — mVt)-!- = — i — u r , (16) 
Bq r 

b v dfl b r du r ( d 



i{u - mn) T = r ^T —fr-\Tz- 2^J u " (17) 

b Tn 

i(u — mQ)— = —i — u z . (18) 
B r 

Finally, the equation of continuity, equation (11), is reduced to 



i{oj — mVL)hi = —c 2 



du r ( d z 

f 7T" - T7TT \U, 



(19) 



dr \dz H 2 , 

Now, we further simplify equations (15) and (17). The last term, —i(m/r)c\(b z /B ), of 
equation (15) can be expressed in terms of u z by using equation (18). The result shows that 
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the term of —i(m/r)c\(b z / Bq) is smaller than the left-hand term, i(u — mQ)u z , by a factor of 
c\/r 2 Q 2 . Considering this, we neglect the last term on the right-hand side of equation (15). 
Next, we consider equation (17). The first term on the right-hand side, r(dQ/dr)(b r /Bo), is 
smaller than the second term, —du r /dr, by a factor of A/r, which can be shown by expressing 
b r in terms of u r by using equation (16). Hence, we neglect the term in the following analyses. 

After introducing the above approximations into equations (15) and (17), we multiply 
i(oj — rriQ) to equation (15) in order to express hi and b v /Bo in equation (15) in terms of u z 
and u T by using equation (19) and (17). Then, after changing independent variables from (r, 
z) to (r, rj), where t] is defined by r] = z/H, we have 



d 2 d (u-mn) 2 -n 2 ± u2 
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This is the basic wave equation to be solved in this paper. 

3. Nearly Vertical Oscillations 

Equation (20) is now solved by approximately decomposing into two equations describing 
behaviors in vertical and radial directions. 

3.1. Nearly Vertical Oscillations 

As mentioned before, we are interested in nearly vertical oscillations (i.e., vertical p- 
mode). The main terms in equation (20) are thus those of the first brackets, and the terms of 
the second brackets are small perturbed quantities. 2 Although the terms of the second brackets 
are small quantities, they are of importance to determine the wave trapping, as is shown in 
subsequent two subsections. 

First, we should notice that the quantity (u — mQ) 2 — Vl\ depends weakly on radius r. 
Hence, in order to consider this weak r-dependence of (u — mQ) 2 — Q 2 ± as a small perturbed 
quantity, the third term in the first brackets of equation (20) is now expressed as 



(id — mtt) 2 — 2 



c 2 



(uj — mil) 2 — jj2 



c 2 



+ 6(r), (21) 



where the subscript c represents the value at capture radius, r c , which is the outer boundary 
of the propagation region of oscillations and will be determined later, and e(r) is a small 
quantities depending on r. The magnitude of e(r) is found from equation (21) when r c and u 
are determined later by an eigen-value problem in the radial direction (see the final subsection 
of this section). If the term of e(r) is transported to terms of small perturbations, in the lowest 
order of approximations, equation (20) is written in the form: 



In nearly vertical oscillations, we have u z ~ h\/c s [equation (15)]. Since the radial component of equation 
of motion [equation (13)] shows that u r ~ (l/QA)/ii, we have u r ~ (H/X)u z , i.e., u r is smaller than u z by a 
factor of H/X. Hence, the terms with the second brackets of equation (20) is smaller than the terms with 
the first brackets by a factor of (H/X) 2 . 
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where the superscript (0) is attached to u z in order to emphasize that it is the quantity of the 
lowest order of approximations. By imposing the boundary condition that does not grow 
exponentially at z = ±oo, we find that the ^-dependence of can be expressed by a Hermite 
polynomial and the term in the large brackets in equation (22) is determined as the eigenvalue 
and found to be n — 1, where n is a positive integer (Okazaki et al. 1987). That is, we have 

u^ = f(r)g^( V ), (23) 

where 



and 



g i0) (v)=Un-i(ri), n = 1,2,3... 
uj - mQ) 2 - Q 2 



n — 1. 



(24) 



(25) 



Here, it is noted that the eigenfunction is taken to be ~H n -i, not T-L n . The reason is that 
in many previous studies hi is adopted as the dependent variable (not u z ) and the ^-dependence 
of hi is taken to be proportional to H n (e.g., Okazaki et al. 1987). The node number of u z 
in the vertical direction is usually smaller than that of hi by unity (e.g., see the equation of 
the ^-component of euation of motion). Considering this, we have adopted H n -i(r]) instead 
of H n (r)) for u z . As shown in equation (23), in the lowest order of approximations, uf\r,rj) 
is expressed in a separable form with respect to r and rj. The r-dependence of is free at 
this stage, which is denoted by f(r) in equation (23). It will be determined later by solving an 
eigen-value problem in the radial direction. 

Equation (25) can be rewritten in the form 



/ oj — mQ\ 2 
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(n-l) + l. (26) 
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In the limit of c\ = 0, this equation is reduced to (u — mQ) 2 = nfl 2 ^, which is the expected result 
from the local dispersion relation in isothermal disks, i.e., [(u — mVt) 2 — h 2 ][(uj — mQ) 2 — nQ 2 ± ] = 
c 2 k 2 (u — mVt) 2 (Okazaki, et al. 1987), where k is the radial wavenumber of oscillations. This 
dispersion relation shows that in non-magnetized isothermal disks, the propagation region of 
the nearly vertical oscillations is described by (to — mQ) 2 > nVt\. This means that for oscillations 
with u, one of their propagation region is uj < mVL — n 1 ^ 2 ^}^. That is, the outer boundary of the 
propagation region on the uj - r plane is given by uj = mVL — n l / 2 Vt^_. As shown later, equation 
(26) suggests that in the present magnetized disks, the outer boundary of the propagation 
region is given by 



uj = mQ- 
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This is really the case as is shown later. 



3.2. Derivation of Equation Describing Radial Behavior 

Now, we proceed to take into account the deviation of the oscillations from purely vertical 
ones as perturbations. We see soon that separation of u z into two functions of r and rj is no 
longer valid. Hence, we consider the effects of small perturbed quantities by introducing a weak 
r-dependence in g. That is, u z is now written as 

u z (r, V ) = f(r)[g {0) (v) + + •••]■ (28) 

Then, from equation (20) we obtain, as an equation describing fg^ 1 ', 
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where is the lowest order expression for u r (r,r/). 

The next subject is to solve equation (29). To do so, 
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In the lowest order of approximations of nearly vertical oscillations, the term of du r /dr on the 
right-hand side of equation (19) can be neglected in evaluating hi, compared with the term of 
(d/dz — z/H 2 )u z . Hence, by using equation (19) we can express dh^ /dr on the right-hand 
side of equation (30) directly by u z . Furthermore, the main term on the right-hand side of 
equation (17) is the term with u^ z \ Hence, by using equation (17) the term of d(b^ / B) / 'dr on 
the right-hand side of equation (30) can be also expressed in terms of . Consequently, in 
the lowest order approximations of nearly vertical oscillations, can be expressed in terms 
of alone from equation (30). After some manipulations we have finally 

d 



u 



(0) 



.drj 



where C s and C\ are operators defined by 



and 



C: 



cl/H 



(oo — mVt) 2 + k 2 



cl/H 



-(tjj-mVL) 2 + K 1 



d_ 

dr 

d_ 

dr 



dr 



d\n(u — mQ) 

dr 



(31) 



(32) 



(33) 



Now, we return to equation (29). The equation is an inhomogeneous equation with 
respect to g^(r,rj). The right-hand side of the equation is now expressed in terms of fg^ by 
using equation (31). As is done in a standard perturbation method, g^(r,rj) is now expressed 
in a series of orthogonal functions of the zeroth order equations as 

9 {1) (r,r ] ) = J2Cm(r)n m (v)- (34) 
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The quantity e(r) is then obtained from the solvability condition of equation (29), using the 
orthogonality of the Hermite polynomials, which is 



+ H^M)(n^ - * (j- - \,)H^) = 0, (35) 

where (A{ji)B{rj)) is the integration of A{rf)B{7]) with respect to rj in the range of (— oo, oo) 
with the weight exp(— t] 2 /2). 

This solvability condition leads to an ordinary differential equation of f{r), when the 
integration with respect to 77 is performed. After some manipulation we can write the results 
in the form: 



K 2 + ^cfl 



l 2 \ d 



oj — mfl d ( f(r 



+ e(r)/(r) = 0. (36) 



. (lo — mil) 2 — k 2 dr \u — mVt,, 
In the previous studies (paper I) on nearly vertical oscillations in non-magnetized disks, we 
have adopted h\ (not u z ) as a dependent variable. To compare our present results with those 
in the previous ones, we introduce here a new variable / defined by / = / /{u — mf2). 3 Then, 
equation (36) is reduced to 



d 



oj — mQ dr 



uj — mVt df 



(u — mil) 2 — k 2 dr 



nn\H 2 



7 = 0. (37) 



In the limit of c\ = 0, this equation becomes formally the same as that used in paper I. 4 

3.3. Radial Eigenvalue Problems 

Next, we solve equation (37) as an eigen- value problem to know where the oscillations 
are trapped and how much the eigen-frequency of the trapped oscillations are. The same WKB 
procudures as Silbergleit et al. (2001) used are adopted here (see also paper I). That is, we 
introduce a new independent variable r(r) defined by 
r r Cj 2 (r') - K 2 (r') 



t r 



-dr', r c = r(r c ), (38) 



where u is defined by Co = (u — mQ), and r\ is the inner edge of disks where a boundary condition 
is imposed. Then, equation (37) is written in the form: 



In the lowest order of approximations, the equation of continuity gives — m£i)h\ + (c^/ H)(d/drj — r])u z = 0. 
If u z is taken to be proportional to H„-i(?7), i.e., u z = /(r)H„_i(?7), the above continuity relation shows 
that h\ has a component proportional to T-L n (rr), i.e., h\ = fh^T-iniv) and f(r) and fh{r) is related by 
i(Lj-mQ)f h = (4/H)f. 

In paper I, polytropic disks are considered. Hence, even in the limit of c\, equation (37) does not become 
idcntial with equation (33) in paper. 
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where 



+ O/=O, (39) 
« T) = ^OT (40) 



Equations (39) and (40) show that the propagation region of oscillations is the region 
where Q > 0, which is the region of e > 0. The region of e(r) > is found to be inside of r c from 
the following considerations. Let us assume tentatively that r c is at a certain radius, although 
it should be determined after solving equation (39). Since we take so that equation (25) holds 
at r c , equation (21) defining e gives 

J i ^2 -| 1/2 

n ± . (4i) 



oj = mVt- 



n-l) + l + e 



cl + cl/2 



This equation gives the u - r relation for a given e. The u - r relation for e = is equation (27) 
and is shown in figures 1 and 2. The curve monotonically increases inwards as r decreases (see 
figures 1 and 2). This means that if r decreases from r c keeping e = 0, the frequency given by 
equation (27) becomes larger than the frequency determined by r c . Hence, to satisfy equation 
(41) under keeping u at the value determined by r c , we must take a positive e. If r increases 
from r c keeping u, on the other hand, equation (41) can be satisfied by taking a negative e. In 
summary, we have e > inside of r c , while e < outside of r c . 5 

We solved equation (39) by a standard WKB method with relevant boundary conditions 
(for details, see Silbergleit et al. 2001). The WKB approximation shows that the solution of 
equation (39) can be represented as 

/ oc g~ 1/4 (r)cos [$(r) - $ c ] (42) 

in the whole capture region < r < r c , except small vicinities of its boundaries of r = (i.e., 
r = r;) and r = r c (i.e., r c ). Here, $(r) is defined by 



$( T ) = f T Q^(r')dr' = f ' Q^{r' f " ^ V , (43) 

Jo Jri —u(r) 

and $ c is a constant to be determined by boundary conditions. To determine the outer boundary 

condition, we take into account the fact that the capture radius, r c , is a turning point of equation 

(39) since the sign of e changes there. The inner boundary condition we adopted is / = or 

df jdr = at r;. As the inner boundary radius we take the marginary stable radius. Then, 

WKB analyses show that the trapping condition is 

i/2 dT = { < n r + V4) for df/dr = 
o V I 7r(n r + 3/4) for 7 = 0, K ' 



It is noted, however, that in the case where the curve of ui = mfl — [{n — l)(c^ + c\)/{c^ + 0.5cjJ + 1] 1 / 2 Oj_ 
monotonically decreases inwards, the situations are changed (this is realized, for example, in the case where 
n = 3 and c\/c^ > 2). That is, we have e > outside of r c and e < inside of r c . This means that a 
propagation region of waves is outside of r c and the waves are not trapped. 
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where n r (= 0,1,2,...) is zero or a positive integer specifying the node number of / in the radial 
direction. The constant <3> c is determined as 



For a given set of parameters, including spin parameter a* and mass of neutron stars, 
M, any solution of equation (44) specifies r c , which gives u of the trapped oscillation through 
equation (26). In other words, u and r c are related by equation (26), i.e., u = u(r c ) or r c = r c (u), 
and the trapping condition determines r c or to as functions of such parameters and 



4. Numerical Results 

To obtain numerical values of the frequency, u, and the capture radius, r c , of trapped 
oscillations, we must specify the radial distribution of acoustic speed, i.e., c s0 (r). The final 
results of numerical calculations show that the trapped region is in the inner region of the 
disks. Hence, we consider the temperature distribution in the standard disk where gas pressure 
dominates over radiation pressure and opacity mainly comes from the free-free processes, and 
adopt (e.g., Kato et al. 2008) 



where is the Eddington luminosity. Parameters a and m affect on the frequencies of trapped 
oscillations only through the magnitude of c s0 . We adopt, throughout this paper, a = 0.1 and 
rh = 0.3. A parameter specifying the strength of magnetic field is c 2 A /c 2 . In this paper, we 
consider the disks where the parameter c\/ c 2 is in the range of c\/ c 2 = ~ 2. Other parameters 
specifying the disk-star system are m(= M/M & ) and a*. We consider the cases of M/M & = 2.0 
and a* = ~ 0.3. 

We only consider two- armed oscillations with one, two, or three node(s) in the vertical 
direction, i.e., n = 1, 2, or 3. 7 Oscillations with more nodes in the vertical direction are less 
interesting from the view point of observability. The inner boundary of oscillations is taken at 
the radius of k = 0, i.e., at the radius of the marginally stable circular orbit. In this paper, 

6 In this section and hereafter, m is often used to denote M/Mq without confusion with the azimuthal 
wavenumber m of oscillations. 

7 In oscillations with n, we have h\(r, z) oc T-L n {z/ H) and u z (r,z) oc TL n ^i(z / H). That is, in oscillations with 
n = 3, u z is plane-symmetric with respect to the equatorial plane, and has one node (where u z — 0) above 
and below the equator. 




for df jdr = 
for / = 0. 



(45) 



M. 
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u z = (i.e., / = 0) is adopted at the radius as a boundary condition, except in figure 5. In 
figure 5, boundary condition of du z /dr (i.e., df jdr ~ 0) is considered as well as u z = in order 
to see effects of boundary condition on results. We find that the differences of the boundary 
condition bring about quantitative differeces in results, but there is no essential differences in 
parameter dependences of results. Hence, except in figure 5, we adopt u z = as the inner 
boundary condition throughout this paper. The horizontal node number, n T , of oscillations we 
consider is mainly n r = and supplementally n r — 1 and 2. 

Figures 1 and 2 are the propagation diagrams for oscillations of n = 1 and 2 (figure 1) 
and n — 3 (figure 2), respectively, in the disks with c\/c 2 = 1 and a* = 0. Only the oscillations 
of n r = are shown in figure 1, but three modes of oscillations, i.e., n r = 0, 1, and 2, are shown 
in figure 2. The propagation regions of oscillations on the frequency-radius diagram is below 
the curve given by uj = mQ — [(n — l)(c 2 + c\)/(c 2 s + c\/2) + l] 1 / 2 !^ [see equation (27)]. The 
results of numerical calculations show that the oscillations of n = 1 with n r = are trapped in 
the radial range shown by the upper thick horizontal line in figure 1. The frequency uj and the 
capture radius r c are, respectively, uj = 798Hz and r c = 3.72r g . Outside r c , the oscillation is 
spatially damped. The radial range of trapped oscillations of n = 2 with n T = is shown by the 
lower thick horizontal line in figure 1. The frequency and the capture radius in this case are 
uj = 359Hz and r c = 3.84r g . 

Trapped oscillations of n = 3 have frequencies lower than those of n = 1 and 2, since on 
the propagation diagram the curve of u — mVL — [(n — l)(c 2 + c^)/(c 2 + c\/2) + l] 1 / 2 ^ is below 
those in the cases of n = 1 or n = 2. In figure 2, the frequency and the radial extend of trapped 
oscillations of n = 3 are shown for three modes in the radial direction; the fundamental mode 
(i.e., n r = 0) and the first two overtones (i.e., n v = 1 and 2). The sets of frequency and capture 
radius for these three modes of n r = 0, 1, and 2 are, respectively, (49.5Hz, 4.60r g ), (28.4Hz, 
6.66r g ), and (17.2Hz, 9.31r g ) in the disks with c^/c 2 = 1.0 and a* = 0. 

Figure 3 shows the c^/c 2 - dependence of the capture radius r c . As a typical case, the 
dependence is shown for oscillations of n v = 0. No spin of the central source is adopted. It is 
noted that when c^/c 2 is close to 2, the capture radius of oscillations of n = 3 is far outside 
and the frequencies are low. These characteristics become more prominent for oscillations with 
n r > 1, although they are not shown in figure 3 (see figure 5). As c^/c 2 increases beyond 2, 
the oscillations of n = 3 are no longer trapped. This is related to the behavior of the curve of 
uj = mVL - [{n- l)(c 2 + c|)/(c 2 + c\/2) + l] 1/2 fi ± on the uj - r plane. If n = 3 and c\/c 2 s > 2, uj 
given by the above relation is negative and is a monotonically increasing function outwards on 
the uj - r plane. Then, the region of e > (i.e., propagation region) is in the outer region and 
there is no trapping (see the previous section). 

The frequency- c^/c 2 relations are summarized in figure 4 for two disks with a* = and 
a* = 0.2. Modes of oscillations adopted are n = 1, 2, and 3. In all cases n r is taken to be n T = 0. 
As mentioned before, the oscillations with n = 3 have low frequencies. In order to examine 
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Fig. 1. Frequency - radius plane (i.e., propagation diagram) showing the propagation region of two-armed 
(to = 2) nearly vertical oscillations (i.e., vertical p-mode oscillations) in vertically isothermal disks with 
toroidal magnetic fields. The value of c\/c 2 = 1.0 has been adopted. The propagation region of oscillation 
modes with n is below the boundary curve labelled by mfl - [(n - 1 ) (eg + c% ) / (c s 2 + . 5c\ ) + 1] 1 / 2 fl ± . The 
boundary curve is shown for two cases of n = 1 and n = 2. In the case of n = 1 (and n r = 0), the capture 
(trapped) zone and the frequency of the trapped oscillations are shown by the upper thick horizontal line 
(the frequency uj is ~ 798Hz and capture radius r c is ~ 3.72r g ). In oscillations with n = 2, the trapped 
oscillation with n r = is shown by the lower thick horizontal line. The trapped frequency uj is ~ 360Hz 
and r c ~ 3.84r g . The inner boundary condition adopted at n is / = 0. This inner boundary condition is 
adopted in all cases in this paper, except for in figure 5. The central star is assumed to have no spin. The 
mass of the central star is taken to be 2Mq in all cases shown in figures in this paper. 

characteristics of these low frequency oscillations more in detail, the frequency-c^/Cg relation 
in case of n = 3 is again shown in figure 5, including cases where other parameter values are 
adopted. That is, in addition to oscillations with n r = 0, oscillations with n r = 1 and 2 are 
considered in figure 5. In addition, the cases where (du z )/dr = is adopted at r ; as the inner 
boundary condition are shown by thin curves. In figure 6, the frequency - spin relation is shown 
for three modes of oscillations with n — 1,2, and 3, where n r = and c\/c^ = 1.0 are adopted. 

5. Discussion 

In this paper we have examined characteristics of trapping of two-armed (m = 2), nearly 
vertical oscillations (vertical p-mode oscillations), assuming that the disk is isothermal in the 
vertical direction and is subject to purely toroidal magnetic fields. Fot mathematical simplicity, 
the ratio of the Alfven speed to the acoustic apeed is constant in the vertical direction, i.e., 
c^/Cg is constant in the vertical direction. The effects of magnetic fields on characteristics of 
the oscillation are measured by taking c\/c^ as a parameter. 

A differece from paper I is that in paper I we have considered a polytropic gas (not 
an isothermal gas) and examined how the characteristics of trapping of the nearly vertical 
oscillations depend on the polytropic index characterizing the gas. In this paper, the disk gas 
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Fig. 2. The same as figure 1, except that the oscillations with n = 3 are considered here. The propagation 
region of the oscillations is below the curve labelled by mfi - [(n - 1) (c 2 + c 2 A )/(c 2 + 0.5c\) + 1] ^ 2 fl x , where 
n = 3 is taken. Trapping of three modes of oscillations with n T — 0, 1, and 2 are shown by three horizontal 
thick lines. The sets of frequency and capture radius for these three oscillation modes are, in turn, (49.5Hz, 
4.60r g ), (28.4Hz, 6.66r g ), and (17.2Hz, 9.30r g ). It is noted that c\/c 2 = 1.0 is adopted here, but in the 
case of c^/Cg = 2.0, there is no trapped oscillations. 




CaVCs 2 

Fig. 3. Capture radius, r c , as functions of c\/c 2 for three modes of n = 1, 2, and 3. The radial node 
number, n r , is taken to be zero with boundary condition / = at n = 3r g , the spin parameter being 
a* = 0. In oscillations with n = 3, the trapping zone extends infinity as c\/ 'c 2 approaches c\/c 2 = 2, and 
no trapping for c\/c 2 > 2. No spin of the central source is taken, i.e., a* = 0. 
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Fig. 4. Frequency-c^/Cg relation of trapped oscillations for some values of vertical node number n and 
spin parameter a*. The oscillations with no node in the radial direction (n r = 0) are considered. 



250 




c A 2 /c s 2 

Fig. 5. Frequency- c^/Cg relation for oscillation modes with n = 3. This figure is to demonstrate the 
effects of differences of radial node number n r and of boundary condition on the frequency of trapped 
oscillations. Two cases of boundary conditions, f\ = [which corresponds to (ft-i)i = 0] and (df/dr)i = 
[which corresponds to {dh\/dr)i =0] are compared for three modes of oscillations with radial node number 
n r — 0, 1, and 2. The thick curves are for the cases where the inner boundary condition is taken as /; = 0, 
while the thin curves are for the cases of (df /dr)\ = 0. The spin parameter a* is taken to be zero. 
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Fig. 6. Frequency-spin relation for modes with n=l, 2, and 3. The radial node number is taken to be 
n r = 0, c^/cg = 1.0, and a* = have been adopted. 

is taken to be isothermal, for simplicity, differet from the paper I, but is subject to toroidal 
magnetic fields. Under these situations, we examine how the characteristics of nearly vertical 
oscillations change by a change of strength of magnetic fields. 

In this paper a partial differential equation [equation (20)] has been solved by a per- 
turbation method. In this procedure, we have assumed that e(r) defined by equation (21) is 
a small positive quantitiy (i.e., e < 1) in the wave propagation region. The final results show 
that this is really acceptable as a first step to examine qualitative behavior of trappeing. To 
do quantitative argument, however, the approximation should be improved especially in the 
case of n = 1 and n = 2. That is, e(r) is zero at r = r c by definition and increases inwards 
monotonically and becomes a maximum at r = r\. The results of calculations show that the 
maximum value of e is 0.46 for n = 1, 0.35 for n = 2, and 0.12 for n = 3 when c\/c^ = 1.0, a* = 
and n r = 0. The value slightly increases with decrease of c\/c^ and increase of a* and n r . 

In this paper we did not quantitatively consider the effects of c s o(r) on frequency. An 
increase of c s0 without any change of other parameters leads to decrease of frequency of trapped 
oscillations. The reason is that an increase of c s o decreases Q. Hence, to satisfy the trapping 
condition (44), an increase of r c is necessary, which leads to decrease of frequency (see figures 
1 and 2). The effects of changes of other various parameter values on frequencies of trapped 
oscillations are qualitatively the same as those in paper I (see tablel in paper I). 

The purpose of this paper is to demonstrate the importance of two-armed (m = 2), 
nearly vertical oscillations as one of possible candidates of disk oscillations describing quasi- 
periodic oscillations observed in low- mass X-ray binaries (LMXBs). One of reasons why we 
take attention on these oscillations is that they can be trapped in the inner region of disks 
and their frequencies can cover a wide range of frequency by (i) differece of modes (n and n r ), 
and (ii) change of disk structure. As a change of disk structure, we considered a change of 
polytropic index in paper I, while we consider here a change of magnetic fields. 
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The parameters describing the difference of oscillation modes are n{= 1,2,3,...) and 
n r {= 0,1,2,...), where n and n r are, respectively, the node number of hi in the vertical and 
horizontal directions. As in paper I, the trapped oscillations of n = 1 and n = 2 have frequencies 
on the order of kHz QPOs, while those of n = 3 have lower frequencies and on the order of 
horizontal branch and normal branch oscillations. An interesting result obtained in this paper is 
that the frequency of trapped oscillations decreases as magnetic fields become stronger (figures 
4 and 5). This is particularly so for oscillations of n = 3 (figure 5). That is, the trapped region 
of oscillations of n = 3 becomes wide and their frequencies decreases as c\/ c 2 increases, and 
finally there is no trapping for c\/c 2 > 2 (figure 5). This situation is similar to the case where 
the polytropic index is increased in polytropic gases (see Paper I). 

The decrease of frequency of trapped oscillations with increase of c\/c 2 can be under- 
stood by considering the following situations. Since toroidal magnetic fields are considered here, 
the nearly vertical oscillations are the fast mode among three MHD oscillation modes. Their 
frequencies seen from the corotating frame are thus higher than those of pure acoustic oscilla- 
tions in non- magnetized disks. If we consider purely vertical oscillations, neglecting u r and u^, 
in vertically isothermal disks with c\/c 2 = const., we find easily that the eigen-frequency of the 
oscillations in the corotating frame, (oj — mf2), is given by 



[see equation (27)], the right-hand side being larger than nVt\ by the presence of toroidal 
magnetic fields, as expected. If the terms neglected in deriving equation (48) are taken into 
account, (u — mQ) 2 becomes larger than the right-hand side of equation (48). Hence, we see 
that one of the propagation region of the nearly vertical oscillation is given by 



This shows that the upper boundary curve specifying the propagation region of oscillations, 
i.e., u = mVL — [(n — l)(c 2 + c\)/(c 2 + 0.5c\) + l] 1//2 f2j_, moves downwards on the propagation 
diagram (u - r plane) as c\/c 2 increases (see figures 1 and 2). In the case of the oscillations 
of n = 3, the boundary curve tends to close to oj ~ (as c\/c 2 approaches 2.0), since the term 
inside the brackets of equation (49) becoms 4. This is the reason why we have no trapped 
oscillations for c\/c 2 ~ 2 when n = 3 and a* = 0. 

Frequencies of QPOs observed in LMXBs have time change, distinct from those in galac- 
tic black hole candidates. So far as trapped oscillations are concerned, a large change of their 
frequencies cannot be expected, unless magnetic fields are considered. Global toroidal magnetic 
fields in accretion disks are time-dependent by amplification due to winding and by damping 
due to reconnection. In the present model of QPOs, the time change of QPOs is attributed to 
the time change of global magnetic fields, but careful discussions are necessary whether this is 
consistent with observations. 




(48) 




(49) 
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There are some important problems remained to be clarified. First, it is not clear 
whether the innermost part of disks can be regarded as a boundary where oscillations are 
reflected back outwards. In disks of standard or ADAF disks, the innermost part of disks will 
reflect incomming waves, since the density decreases there sharply inwards. [See Kato et al. 
(1988) and Manmoto et al. (1996a,b) for reflection of waves in the inner edge of disks.] In slim 
disks, however, there will be no sharp density decrease inwards near the transonic radius, and 
thus the reflection of waves in the innermost region will be not so efficient. In the case where 
the central source is a neutron star, the stellar surface or a transition region near to the surface 
will, at least, partially reflect incoming waves. 

Whether nearly vertical oscillations are really excited on disks is also a problem to be 
examined, since they will not be excited by thermal and viscous overstable processes (e.g., Kato 
1978), because of the presence of node(s) in the vertical direction. Most conceivable processes 
will be stochastic processes of turbulence (Goldreich and Keely 1977a,b). In many stars with 
various characteristics (e.g., different evolutionary stages, effective temperatures...), solar-like 
(non-radial) oscillations have been observed. The so-called K-mechanism cannot excite all of 
these oscillations. Their origin is now known to be stockastic processes of turbulence [see 
Samadi (2009) for a recent review of stochastic excitation of the oscillations]. Compared with 
in stars, much stronger MHD turbulence by magneto-rotational instability (MRI) are expected 
in accretion disks, especially in the inner region of disks. Hence, it is natural to suppose that in 
accretion disks many trapped oscillations are simultaneously excited by turbulence. This may 
be one of causes of vareity of QPOs in LMXBs. This is a problem to be examined in the future. 
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